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Abstract 

We identify within the SU (5) framework the minimum number of soft supersym- 
metry breaking parameters which can yield a bottom squrak (sbottom) as the next to 
lightest supersymmetric particle. We focus in particular on the neutralino-sbottom 
coannihilation scenario which gives rise to the desired neutralino dark matter relic 
density. We find solutions in which the sbottom mass is greater than or of order 
210 GeV, while the gluino and the first two family squarks are heavier than 1 TeV. 
Some benchmark points which can be tested at the LHC are presented. 
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1 Introduction 

In a recent paper [1], the authors have shown that b-r Yukawa coupling unification, 
realized in some well motivated SU{5) and SO (10) models, is compatible with the 
constrained minimal supersymmetric model (CMSSM) and with neutralino dark mat- 
ter abundance only if there exists neutralino-stop coannihilation. In order for this 
coannihilation scenario to be effective, the lighter stop must be the next to lightest 
supersymmetric particle (NLSP), and it is quasi-degenerate in mass with the dark 
matter lightest supersymmetric particle (LSP) neutralino (to within 10-20% or so.) 
A very recent analysis [2] of this neutralino-stop cooannihilation scenario shows that 
the ATLAS search for supersymmetry [3], corresponding to an integrated luminos- 
ity of 1 fb _1 , essentially rules out an NLSP stop mass below 140 GeV. Future LHC 
searches will no doubt provide far more stringent constraints on this scenario, or 
perhaps discover neutralino-stop coannihilation. Another colored particle which can 
have mass of order 100 GeV and can still survive the recent ATLAS bounds [4] is the 
sbottom quark when it is quasi-degenerate in mass with neutralino. 

Motivated by these considerations we propose to investigate the neutralino-sbottom 
coannihilation scenario in this paper. The analysis in [1] shows that this scenario 
requires a framework larger than the CMSSM with additional soft supersymmetry 
breaking (SSB) parameters. We have selected to investigate neutralino-sbottom coan- 
nihilation in SU(5), which naturally allows for the presence of additional parameters 
in accord with minimal supergravity [5] . We should note here that neutralino-sbottom 
coannihilation in SU(J5) has previously been explored in the references listed in [6]. 
Our analysis, we believe, provides for the first time a comprehensive study of this 
scenario in SU(5). 
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The outline for the rest of the paper is as follows. In section 2 we briefly describe 
the model, list the 577(5) inspired SSB parameters, and the range of values employed 
in our scan. Section 3 describes the scanning procedure and the relevant experimental 
constraints that we have employed. The results pertaining to neutralino-sbottom 
coannihilation are discussed in section 4, and our conclusions are summarized in 
section 5. 

2 Soft Supersymmetry Breaking (SSB) Parame- 
ters in SU(b) 

We will search for neutralino-sbottom coannihilation in the SU (5) framework and, for 
simplicity, we assume that the SSB mass terms for sfermions are family independent. 
In SU (5) the standard model (SM) fermions per family are allocated to the following 
representations: 5 D (d c ,L) and 10 D ((3,^ c ,e c ), where in the brackets, we have 
employed standard notation for the SM fermions. It seems natural to consider two 
independent SSB scalar mass terms, at Mq, and mio, for the matter multiplets. 
The MSSM Higgs doublets belong to 5(H U ) and h(Hd) representations of 577(5), to 
which we assign two independent SSB mass terms, vrtu u and rriH d - 

The minimal SU (5) model predicts b-r Yukawa coupling unification at Mq from 
the interaction 5.10.5# d , and its LHC implications have been discussed in [1]. It was 
shown in [10] that b-r Yukawa coupling unification can be relaxed by including either 
non-renormalizable interactions in the theory, or by employing a more complicated 
Higgs sector, beyond the minimal one. To implement neutralino-sbottom coanni- 
hilation, we find it helpful not to require b-r Yukawa unification. We also invoke 
non-universal soft trilinear terms. In particular, we assume that at Mq, there are 
two independent parameters A^(— A T ) and A t) where the equality = A T keeps the 
number of SSB parameters to a minimum. To summarize, we consider the following 
SSB terms for this study: 

Mi/2, m 10 ,m^m Hu ,m Hd , tan fi,A t ,A h (= A T ) and sign(n). (1) 

We shall see that this is the minimal set of independent MSSM terms which allows 
one to implement neutralino-sbottom coannihilation in SU(5). We will set fi > in 
this paper. 

3 Phenomenological constraints and scanning pro- 
cedure 

We employ the ISAJET 7.80 package [13] to perform random scans over the pa- 
rameters listed in Eq.(l). In this package, the weak scale values of gauge and third 
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generation Yukawa couplings are evolved to Mq via the MSSM renormalization group 
equations (RGEs) in the DR regularization scheme. We do not strictly enforce the 
unification condition g^ — g\ — #2 at Mq, since a few percent deviation from unifica- 
tion can be assigned to unknown GUT-scale threshold corrections [14]. The difference 
between gi(= #2) an d #3 at Mq is no worse than 4%. 

The various boundary conditions are imposed at Mq and all the SSB parameters, 
along with the gauge and Yukawa couplings, are evolved back to the weak scale 
Mz- In the evaluation of Yukawa couplings the SUSY threshold corrections [15] are 
taken into account at the common scale Msusy = -\/ m t L m t R i where mi L and m^ R 
are the third generation left and right handed stop quarks. The entire parameter 
set is iteratively run between M z and Mq using the full 2-loop RGEs until a stable 
solution is obtained. To better account for leading-log corrections, one-loop step-beta 
functions are adopted for gauge and Yukawa couplings, and the SSB parameters 
are extracted from RGEs at multiple scales = m^(m^). The RGE-improved 1- 
loop effective potential is minimized at an optimized scale Msusy, which effectively 
accounts for the leading 2-loop corrections. Full 1-loop radiative corrections are 
incorporated for all sparticle masses. 

The requirement of radiative electroweak symmetry breaking puts an important 
theoretical constraint on the parameter space. Another important constraint comes 
from limits on the cosmological abundance of stable charged particles [16]. This 
excludes regions in the parameter space where charged SUSY particles, such as f\ or 
ii, become the LSP. We keep only those solutions where the lightest neutralino is the 
LSP which, in most cases, saturates the WMAP (Wilkinson Microwave Anisotropy 
Probe) dark matter relic abundance bound. Neutralino-sbottom coannihilation plays 
an important role in realizing the desired LSP relic abundance. 

We have performed random scans for the following parameter range: 

< m b ,m Hd < 5TeV 
0<m 1Q ,m Hu ,< lOTeV 
< Mi/a < 2 TeV 
-15TeV < A t < 15 TeV 
-15TeV < A b = A r < 30 TeV 
1.1 < tan/3 < 60 

fj, > (2) 

where m t = 173.3 ± 1.1 GeV [17] is the top quark pole mass. We use m^mz) = 2.83 
GeV which is hard-coded into ISA JET. The various boundary conditions in Eq.(2) 
are implemented through the scanning procedure. We found, for instance, that 7% 
has to be smaller than mio to realize the neutralino-sbottom coannihilation solutions. 
In scanning the parameter space, we employ the Metropolis-Hastings algorithm 
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as described in [18]. All of the collected data points satisfy the requirement of ra- 
diative electroweak symmetry breaking (REWSB), with the lightest neutralino in 
each case being the LSP. After collecting the data, we impose the mass bounds on 
all the particles [16] and use the IsaTools package [19] to implement the following 
phenomenological constraints on points that have sbottom coannihilation solution: 

m h (lightest Higgs mass) > 114.4 GeV [20] 

BR(B S < 1.1 x 1(T 8 [21] 

2.85 x 1(T 4 < BR(b -> sj) < 4.24 x 10~ 4 (2a) [22] 

0.15 < gg^g^MssM < 2A1 / 3a A 23 

tt CDM h 2 = O.llllffi (5^) [24] 

As far as the muon anomalous magnetic moment is concerned, we only require that 
the model does no worse than the SM. 



4 Sbottom Coannihilation 

To realize the desired LSP dark matter relic abundance through neutralino-sbottom 
coannihilation, the mass difference between these two should be less than or of order 
20% [6]. In Figure 1 we present our results in the M\/2 - mio/m^ At/A^ - m^jm^ 
At/A}, - nriH u /mH d) m§/rriH d - ^HuI^Ha planes, with the parameter values all defined 
at Mq. Gray points are consistent with REWSB and LSP neutralino. The orange 
points satisfy, in addition, the particle mass bounds, constraints from BR(B S — » 
/i + //~), BR(B U —> tv t ) and BR(b —> sj). The blue points form a subset of orange 
points and correspond to NLSP sbottom, but which is not closely degenerate in 
mass with the neutralino. The red points form a subset of blue points and represent 
neutralino-sbottom coannihilation, the scenario we are really after ! The vertical 
and the horizontal dashed lines show equality of parameters along the x and y-axes 
respectively. 

In order to explain our findings we consider the one loop renormalization group 
equations for the third generation squarks and sleptons [25] 

16712 Jt m l s = ^ + ^-f^W-6^ 2 |M 2 | 2 -^ 1 2 |M 1 | 2 + ^, (3) 

l^ 2 f t m% = 2X t -fg 2 3 \M 3 \ 2 - 3 ^g 2 \M 1 \ 2 -^g 2 1 S, (4) 



167r 2 -m| = 2X b --g 2 3 \M 3 \ 2 --g 2 l \M 1 \ 2 + -g 2 1 S, (5) 
d o „ d i -a i o 6 Oi . r in 3 



167r 2 -m| 3 = X T -6g 2 \M 2 \ 2 --g 2 \M 1 \ 2 --g 2 S, (6) 
167r 2 |m| = 2X T -fg 2 \M 1 \ 2 + 6 -g 2 S. (7) 
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Here 



X t = 2\y t \\m] Iu +ml 3 +ml l ) + 2\A t \\ (8) 

X h = 2\y h \\m] {d +ml 3 +m\^ + 2\A h \\ (9) 

X T = 2\y T \\m 2 Hd +m\ 3 +ml l ) + 2\A T \\ (10) 

S = Tr[l>ijj = - m^ d + (11) 

and Q3, S3, 6/3, L 3 , denote the third generation squarks and sleptons. Also, gi and 
Mi (i = 1,2,3) denote the gauge couplings and gaugino masses for [7(l)y, SU(2)l 
and SU(3)c, and y^, (j = i, 6, r) are the third family Yukawa couplings and 
trilinear scalar SSB couplings, respectively. 

In the M1/2 - ^10/^5 plane of Figure 1, we show that in order to have sbottom 
NLSP, the ratio 777,10/777,5 should be greater than 1.4 (blue points). If we require 
successful neutralino-sbottom coannihilation, the ratio 777,10/777,5 > 2 (red points). In 
general, there is a two step process for realizing NLSP sbottom. First we consider the 
RGE effects, and in the second step we search for cancellation in the 2x2 sbottom 
mass matrix between the diagonal and off diagonal (—777,5(^6 + /itan/3)) entries. 

Let us see why it is not possible to have sbottom NLSP starting from universal 
sfermion masses (777,5 = m io)- First, we note, that it is difficult to make the sbot- 
tom lighter than the stop. Eqs.(3)-(11) show that the gluino loop contribution raises 
the squark masses (we ignore the hypercharge contribution because it is sufficiently 
small), while the Yukawa and SSB trilinear couplings tend to lower them. The left 
handed squarks are heavier than the right handed ones due to the SU(2) L contri- 
bution. From the QCD point of view the sbottom and stop masses 2 renormalize 
identically. To split them, the loop corrections involving Yukawa couplings must play 
a role. Since the top Yukawa coupling, for most tan/3 values, is larger than the bot- 
tom Yukawa coupling, it is hard to make m~ c < m 2 ~ c through RG running. If the soft 
trilinear couplings A t and A^ are independent of each other, one could make one of 
the sbottom mass eigenvalues lighter than the stop. But in this case we need to make 
sure that the eigenvalues of the stau mass 2 matrix remain positive, after suitable 
cancellation in the sbottom mass 2 matrix is accomplished. It turns out, as we show 
in At/At - raio/mg plane (Figure 1), that 777,10 = 777,5 an d A t = A^ are not compatible 
with neutralino-sbottom coannihilation. It is possible to achieve neutralino-sbottom 
coannihilation with A t = A^ = A T) provided 777,10/777,5 > 2. The ratio 777,10/777,5 takes 
its minimal value around 1.4, which corresponds to A t /Ab < —0.5. Note that for 
M1/2 > ^0 where 777,0 is the universal SSB mass term for sfermions at Mq, it is 
difficult to make the sbottom lighter than the right handed stau because of QCD 
corrections from the gluino loop. Combining these two observations, it is clear that 
we need 777,10 > ^5 for any given value of Mi/ 2 - 

In the Af/Afj - mH u /rriH d plane in Figure 1 we see two viable neutralino-sbottom 
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Figure 1: Plots in the M 1/2 - m 10 /m^ A t /A b - m 10 /m^ A t /A b - m Hu /m Hd) m^/m Hd 
- ^Hu/^Hd planes. Gray points are consistent with REWSB, neutralino LSP. 
The orange points satisfy, in addition, the particle mass bounds, constraints from 
BR(B S —> BR(B U — » rv T ) and BR(b — )■ 57). The blue points form a subset 

of orange points that shows sbottom quark as an NLSP. The red points represent 
a subset of blue points that corresponds to the neutralino-sbottom coannihilation 
scenario. The vertical and horizontal dashed lines show equality of parameters along 
x and y-axes respectively. 
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coannihilation regions which have A t f A b = 1 or mH u /mH d = 1. However, we also see 
that it is not possible to have such solutions with both these ratios simultaneously 
equal to unity. 

In minimal SU (5) the MSSM down type Higgs field resides in the 5 representation, 
and so it is interesting to see whether or not neutralino-sbottom coannihilation or 
sbottom NLSP is possible by setting = rriH d . In the m§/rriH d - mH u /^H d plane 
with 7% = rriH d) neutralino-sbottom coannihilation requires m^/m^ > 8 (red 
points), while sbottom NLSP only needs m^/m^ > 2 (blue points). 

In Figure 2 we present results in the 7% - M 1 / 2: f% - ^io, A t - rriH d - mH u , 
tan/3 - 7% and R^ - tan/3 planes which display the values of the SSB parameters 
required to yield sbottom NLSP/ coannihilation. Gray points are consistent with 
REWSB and neutralino LSP. The orange points satisfy, in addition, the particle 
mass bounds, constraints from BR(B S — > /i + /x~), BR(B U — >> tis t ) and BR{b —> sj). 
The blue points form a subset of orange points and correspond to sbottom NLSP. The 
red points represent a subset of the blue points and represent the neutralino-sbottom 
coannihilation scenario. 

As we can see from the m§ - M 1 / 2 panel, neutralino-sbottom coannihilation prefers 
relatively heavy gauginos (Mi/ 2 > 400 GeV), while can be as light as 500 GeV. 
For 7% lighter than 500 GeV, it is difficult to separate the right handed sbottom mass 
from the left handed tau slepton. In the 7% - mi plane, we see that mi should be 
larger than 2 TeV in order to have neutralino-sbottom coannihilation. 

As noted above, the RG running is not sufficient to realize an NLSP sbottom, 
and some cancellation is required in addition. Consider the plot in A t -Ab plane at 
SUSY scale where we see that neutralino-sbottom coannihilation solutions mostly 
require \Af,\ > \A t \, where \Ab\ can be O(10 — 15)TeV. The off-diagonal entries 
m^Ab + /itan /3) for the sbottom quark mass matrix are of comparable magnitude to 
the diagonal entries and, as a result, we can realize neutralino-sbottom coannihilation. 

The plot in rriH d -mH u (at Mq) plane shows that neutralino-sbottom coannihilation 
requires rriH u > ^H d - This can be understood by considering the one-loop RGEs for 
m 2 Hu and m 2 Hd : 

3X t - 6g 2 2 \M 2 \ 2 - ^IMxI 2 + ^5, (12) 
3X b + X T - 6g 2 2 \M 2 \ 2 - ^IMxI 2 - ^g 2 S. (13) 

Combining with Eqs. (8)-(10) we see that the RG evolution of m 2 Hu and m 2 H depend 
on y t) 7/5, \A}\ and \A t \. Our conclusion from Figure 1 is that neutralino-sbottom 
coannihilation solutions require in most cases \A t \ < which means that Tn 2 Hd 

renormalizes (decreases) more than m 2 Hu . In order to have —Tn'j Id > —m 2 Hu at low 
scale, we need to start in most cases with m 2 H < m 2 H at Mq. 



167T —m H 

dt Hu 

i r 2 ^ 2 
lOTT —m 



dt 



H d 
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Figure 2: Plots in the mi/ 2 - mi - m^, A t - rriH d - m H u , tan (3 - and - 
tan/3 planes. Color coding same as in Figure 1. Lines in black have unit slopes with 
x and y coordinates equal 
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We note from - tan /3 plane that sbottom coannihilation prefers small to mod- 
erately large values for tan /3 (7 < tan/3 < 40). In this region the top Yukawa coupling 
reaches its minimum value at Mq. 

In order to quantify b-r Yukawa unification, we define the quantity R^ as 

R bT = maX ^ } . (14) 

The i?5 r -tan/3 plot in Figure 2 shows that it is difficult to reconcile b-r Yukawa 
unification with the neutralino-sbottom coannihilation scenario. For the parameter 
range given in Eq.(2), b-r Yukawa unification is at the level of 20% or so, which agrees 
with the results in ref. [10]. 

In Figure 3, we show plots in the m~ g - mj c , rrtg - vn~ g - , m^-mg , m^-m^, 
and m^-m^o, planes. The color coding is the same as in Figure 1, except that the 
green points in the - m^o plane form a subset of orange points and satisfy all 
current experimental bounds [3]. 

From the results recently presented by the ATLAS [3] and CMS [7] collaborations, 
our findings are just beginning to be tested. In the mj c - m~ g and - rrtg planes we 
show graphs for the first generation squarks versus gluino mass. From the m d ~ c - rrtg 
plane one sees that the minimum value of mj c in the neutralino-sbottom coannihi- 
lation scenario is around 800 GeV, corresponding to a gluino mass of 1.2 TeV. From 
the plot in - vn~ g plane, we see that is heavy (> 2.2 TeV), and lies well above 
the current experimental bound. 

The plot in m~ hi -m~ g plane shows that the neutralino-sbottom coannihilation solu- 
tions (red points) allow a relatively light sbottom quark, with mass as low as 210 GeV. 
The corresponding gluino mass is around 1.2 TeV. For gluino mass around 3 TeV, 
the corresponding value of mg in this scenario is around 700 GeV. The plot in mg 
- mi plane shows that in neutralino-sbottom coannihilation solutions, the minimum 
(maximum) value of the lighter stop quark mass is of order 300 GeV ( 1.4 TeV.) 

In Figure 3 we also show plots in m~ hi - and m~ hi - m^o planes. In the m~ hi - m^o 
plane, the green points represent solutions that satisfy all the constraint mentioned 
above. Points within 20% of the line with unit slope represent neutralino-sbottom 
coannihilation, while points in blue color represent just NLSP sbottom solutions 
(without coannihilation.) From the graph we see that the minimum value of neu- 
tralino mass in either case close to 210 GeV. The minimum value of stau mass 
corresponding to the coannihilation scenario is about 500 GeV. 

In Figure 4 we show a plot in the mg - BR(B S —> /x + /x~) plane, with the color 
coding the same as in Fig. 1. The vertical dashed line on the left marks the SM 
prediction for BR(B S —> /x + /x~), while the dashed line on the right shows the current 
LHCd+CMS limit. As previously mentioned in section 3, we have applied the latest 
combined LHCd+CMS limit for BR{B S -> /x + /x") of 1.1 x 10" 8 . This can be seen 
as a sharp cutoff in solutions represented by yellow points. We observe that the 
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Figure 3: Plots in the m~ g - mj c , m~ g - m^ c , rrig - mg , m^-m^ , and m^-m^ planes. 
Color coding is the same as in Figure 1 except that green points in m^-m^o plane 
form a subset of orange points which satisfies all current experimental bounds [3] . 
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BR(B S /z+/i~) x 10- 8 

Figure 4: Plot in mi -BR{B S plane. Color coding is the same as in Figure 

1. Left vertical dashed line corresponds to the SM prediction for BR(B S — >> 
and right dashed line shows the current LHCd+CMS limit. 

neutralino-sbottom coannihilation solutions predict a value for BR(B S —> of 
around 4 x 10 -9 , close to the SM prediction of around 3.2 x 10 -9 . This prediction of a 
slight excess in BR(B S for the neutralino-sbottom coannihilation scenario 

should be tested soon. 

In Table 1 we present three characteristic benchmark points which satisfy all con- 
straints presented in Chapter 3. Points 1 and 3 respectively represent the minimum 
and maximum values of the sbottom mass (~ 210 GeV and ~ 820 GeV), correspond- 
ing to neutralino-sbottom coannihilation. Point 2 represents just an NLSP sbottom 
solution (without insisting on correct dark matter relic abundance value), with the 
stop and stau masses relatively close to it. Since the LSP is essentially a pure bino, 
both its spin-independent (~ 10 -13 — 10 _11 pb) and spin-dependent 10 -10 — 10 _7 pb) 
cross sections on nucleons are rather small [26]. Consequently it won't be easy to 
detect the LSP in direct and indirect experiments. However, as shown in [2] and [27], 
the LHC supersymmetry searches can be exploited to probe cross sections of this 
magnitude. 

5 Conclusions 

We have described in detail the conditions under which the neutralino-sbottom coan- 
nihilation scenario can be realized in supersymmetric SU(J5). In particular we have 
identified, for the first time we believe, the minimum number of soft supersymmetry 
breaking parameters that are required in order to have NLSP sbottom and neutralino- 
sbottom coannihilation in SU(5). The coannihilation scenario predicts the existence 
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Table 1: Point 1 shows the minimum value of sbottom mass (213 GeV) for neutralino- 
sbottom coannihilation. Point 2 represents a solution with sbottom NLSP (414 GeV), 
more than 50% heavier than the LSP neutralino. Point 3 represents a neutralino- 
sbottom coannihilation solution with a heavier sbottom mass (823 GeV). 

of relatively heavy (> 1 TeV) gluino and first two generation squark masses. The 
NLSP sbottom in the coannihilation scenario is quasi-degenerate in mass with the 
LSP neutralino, and it can be as light as 210 GeV or so, without running into con- 
flict with the current supersymmetry search at the LHC. We also highlight a few 
benchmark points which can be tested at the LHC. 
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